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1 Introduction 

The fundamental theory of the strong interaction is quantum chromo dynamics 
(QCD), which is a relativistic field theory with local gauge invariance, whose el- 
ementary constituents are colored quarks and gluons. In principle, QCD should 
provide a complete description of nuclear structure and dynamics. Unfortunately, 
QCD predictions at nuclear length scales with the precision of existing (and antic- 
ipated) experimental data are not available, and this state of affairs will probably 
persist for some time. Even if it becomes possible to use QCD to describe nuclei 
directly, this description is likely to be cumbersome and inefficient, since quarks 
cluster into hadrons at low energies. 

How can we simplify this problem to make progress? We will employ a framework 
based on Lorentz-covariant, effective quantum field theory and density functional 
theory. Effective field theory (EFT) embodies basic principles that are common 
to many areas of physics, such as the natural separation of length scales in the 
description of physical phenomena. In EFT, the long-range dynamics is included 
explicitly, while the short-range dynamics is parametrized generically; all of the 
dynamics is constrained by the symmetries of the interaction. When based on a local, 
Lorentz- invariant Lagrangian (density), EFT is the most general way to parametrize 
observables consistent with the principles of quantum mechanics, special relativity, 
unitarity, gauge invariance, cluster decomposition, microscopic causality, and the 
required internal symmetries. 

Density functional theory (DFT), which has been widely used in atomic and 
condensed-matter physics, allows us to describe the nuclear many-body system with 
a universal energy functional that depends on nuclear densities and four-vector cur- 
rents. In principle, knowledge of the full energy functional allows us to calculate 
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any observable for the (zero-temperature) many-body system; moreover, a simpli- 
fied treatment of the functional based on quasi-particle orbitals still provides an 
exact description of bulk nuclear properties and some single-particle observables. 
The great advantage of DFT is that calculations of this subset of properties can be 
made without knowledge of the many-body wave function, or with a simple one. 
Finally, if relevant expansion parameters can be found, the energy functional can be 
truncated to a manageable size, and the accuracy of the truncation can be tested 
quantitatively for the observables in question. 

The basic properties of nuclei provide stringent constraints on any nuclear theory. 
An accurate description of these properties is necessary for any useful predictions 
or extrapolations. We certainly want to reproduce the observed shapes of nuclei: 
the interior density of a heavy nucleus should be relatively constant, there should 
be a well-defined surface, and because of nuclear "saturation", the radius R of a 
nucleus should scale according to R cx I? 1 / 3 , where B = N + Z is the total number 
of neutrons and protons. Moreover, the total energy E of the nucleus should agree 
with the "liquid drop" formula 

E = -a\B + a 2 B 2 / 3 + a 3 Z 2 /B 1 ^ 3 + a 4 (7V - Zf/B + • • • , (1) 

where typical values for the a. L coefficients are given in 0121111]. 

The particle spectrum is determined by the qualitative features of the single- 
particle potential. In nonrelativistic (Schrodinger) terminology, the central potential 
is midway between a harmonic oscillator and a square well; this shape determines 
the ordering of the levels as a function of the orbital angular momentum. (See 
PQ, Figs. 57.1 and 57.2.) In addition, the spin-orbit potential is strong, which is 
instrumental in determining the major shell closures and, hence, the nuclear shell 
model. We will see later how these features are easily reproduced in a description 
based on the Dirac equation. 

These simple nuclear features are the ones we will focus on. We expect that 
they can be adequately described by a single-particle equation with an effective, 
one-body interaction. Such an approach has many names, depending on the system 
being studied and on the practitioner: "shell model" , "mean-field theory" , "Kohn- 
Sham" DFT, etc. Our goal is to correlate (fit) a modest number of nuclear bulk and 
single-particle data and then to predict other, similar data as well as possible. 

1.1 Why Use Hadrons? 

Well, why not? Our focus is on low-energy, long-range nuclear characteristics, and 
all measured observables are colorless. (In fact, most of the observables of interest 
to us are dominated by the isoscalar part of the interaction.) Moreover, hadronic 
variables (baryons and mesons) are efficient, since hadrons are the particles that are 
observed in experiments. Colored quarks and gluons participate only in intermediate 
states, and such "off-shell behavior" is unobservable; by using hadrons, we expend no 
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theoretical effort combining quarks and gluons into color singlets that can actually 
be observed. 

So we pick the most efficient degrees of freedom by choosing hadrons. We will 
have to parametrize the nuclear EFT Lagrangian anyway, since we cannot compute 
its true form from QCD, and hadronic variables, if combined in all forms consistent 
with the underlying symmetries, provide sufficient flexibility for our parametrization. 
We cannot guarantee that a single-particle hadronic approach will be successful in 
describing all of the observables of interest, but we want to see how well we can do. 

1.2 Why Use the Dirac Equation? 

To motivate the Dirac equation as straightforwardly as possible, compare the particle 
spectrum (and fine structure) in a light atom with the spectrum in a heavy nucleus. 
An example of the former is given in while an early example of the latter is 
given in jS], which is reproduced in Fig. 57.3 of The most striking result is that 
it is impossible to draw the atomic fine structure to scale, since the splittings are 
roughly 1/10,000 as large as the major-level splittings (at least for the deeply bound 
atomic levels). In contrast, the nuclear spectrum shows that the "fine" structure is 
really "gross"; the spin-orbit splittings are as large as the major- level splittings to 
within a factor of two! 

The implication is that there must be some relativistic effects that are important 
in nuclei (unlike light atoms), and thus it is much more natural to use the Dirac 
equation to describe the quasi-particle nucleon wave functions. 

1.3 Quantum Hadrodynamics (QHD) 

We will refer to Lorentz-covariant, meson-baryon, effective field theories of the nu- 
clear many-body problem as "quantum hadrodynamics" or QHD [BllffllTllSll^l lTUllTTl 
I12[ll3j. When QHD is applied within the framework of modern EFT and DFT, 
it provides a quantitative description of bulk nuclear properties and the spin-orbit 
force throughout the Periodic Table [9|ll4|lll| . This success arises from the presence 
of large Lorentz scalar and vector mean fields, which imply that there are large rela- 
tivistic interaction effects in nuclei under normal conditions ^2j. There is evidence 
from QCD sum rules that these large fields are dynamical consequences of the under- 
lying chromodynamics |15 | I16|. Moreover, similar relativistic effects are responsible 
for the efficient description of spin observables in medium-energy proton-nucleus 
scattering using the Relativistic Impulse Approximation an d they are consis- 
tent with the major role played by scalar and vector meson exchange in modern 
boson-exchange models of the nucleon-nucleon (NN) interaction [T^]. All of these 
features motivate further investigation into the application of QHD to the nuclear 
many-body problem. 
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2 Effective Field Theory 

A modern discussion of QHD begins by interpreting the Lagrangian as defining a 
nonrenormalizable, effective field theory (EFT) |19|l9j . An effective Lagrangian con- 
sists of known long-range interactions constrained by symmetries and a complete, 
non-redundant set of generic short-range interactions (i.e., "contact" and "gradi- 
ent" terms). The division between "long" and "short" is characterized by the break- 
down scale A of the EFT. While it is not possible at present to derive an effective 
hadronic theory directly from the underlying QCD, the EFT perspective implies 
that this is not necessary. If one constructs a general Lagrangian that respects 
the symmetries of QCD: Lorentz covariance, parity conservation, time-reversal and 
charge-conjugation invariance, (approximate) isospin symmetry, and spontaneously 
broken chiral symmetry, then the EFT is a general parametrization of observables 
below the breakdown scale. 1 

For QHD, we identify A with the mass scale of physics beyond the Goldstone 
bosons (pions); we will see that A ~ 600 MeV. At momenta small compared to A, 
short-distance physics (such as the substructure of nucleons) is only partially re- 
solved and so may be incorporated into the coefficients of field operators organized 
as a derivative expansion. The coefficients of these short-range terms may eventually 
be derived from QCD, but at present, they must be fitted by matching calculated 
and experimental observables. In principle, there are an infinite number of (non- 
renormalizable) terms, but in practice, the Lagrangian or energy functional can be 
truncated to work to a given precision [S]. The EFT is useful if this truncation can 
be made at low enough order that the number of free parameters is not prohibitive. 

The EFT perspective, with the freedom to redefine and transform fields, implies 
that there are infinitely many representations of low-energy QCD physics. But they 
are not all equally efficient or physically transparent. One of the possible choices is 
between Lorentz-covariant and nonrelativistic formulations. Recent developments in 
baryon chiral perturbation theory support the consistency (and utility) of a covariant 
EFT, with Dirac nucleon fields in a Lorentz-invariant, effective Lagrangian density 
[20112 lll2"2"] . A similar framework underlies QHD approaches to nuclei. 

In QHD, the only essential degrees of freedom are the nucleons and pions. Only 
these stable particles can appear on external lines with timelike four-momenta. The 
long-range pion-pion and pion-nucleon interactions are included in a nonlinear real- 
ization of the spontaneously broken SU(2) L x SU(2) fi chiral symmetry, which avoids 
dynamical assumptions inherent in linear representations. These interactions can 
be written down systematically, given an appropriate power-counting scheme, to 
be discussed shortly |S]. Low- mass vector mesons are typically included for phe- 
nomenological reasons, but are not required, since their masses are roughly equal to 
the breakdown scale A; they are absent from point-coupling Lagrangians, for exam- 



1 It is straightforward to include the local U(l) gauge symmetry of the electromagnetic interaction 
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pie [231124] . In descriptions of NN scattering and of nuclear structure and reactions, 
the heavy, non-Goldstone bosons appear only on internal lines (with spacelike four- 
momenta) and allow us to parametrize the medium- and short-range parts of the 
NN interaction, as well as the electromagnetic form factors of the hadrons |19l l5j. 
The heavy bosons are also convenient degrees of freedom for describing nonvanish- 
ing expectation values of bilinear nucleon operators, such as NN and N'j^N, which 
are important in nuclear many-body systems 0013 • This explains why it is useful 
to introduce collective degrees of freedom with other quantum numbers, such as a 
A baryon to incorporate important pion-nucleon interactions |2HI25| . Because one 
must always truncate the EFT Lagrangian, these degrees of freedom can be efficient 
in the many-body problem whether or not they are actually observed as hadronic 
resonances. 

A Lorentz scalar, isoscalar mean field in nuclei is an efficient way to include 
implicitly the effects of pion exchange that are the most important for describing 
bulk nuclear properties. Because the chiral symmetry is realized nonlinear ly, one 
can add to the theory a light scalar, isoscalar, chiral-singlet field with a Yukawa 
coupling to the nucleon, just as in the original Walecka model [Bj. Nonlinear self- 
interactions of this new scalar must be included, with adjustable couplings that 
arise in part from the nucleon substructure. Since the expectation value of the pion 
field in nuclear matter vanishes at the mean-field level, one makes the remarkable 
observation that the mean-field theory (MFT) of the Walecka model is consistent 
with chiral symmetry, provided we think in terms of a nonlinear realization of the 
symmetry. The light scalar, isoscalar field, which is not the chiral partner of the 
pion, plays the same role in the EFT as it does in the Walecka model: it simulates 
important tttt and NN interactions that must be included from the outset to generate 
a realistic description of nuclear matter and nuclei. 

To make systematic calculations, the EFT approach exploits the separation of 
scales in physical systems, with the ratios of scales providing expansion parameters. 
A connection between appropriate QCD scales and nuclear phenomenology is made 
by applying Georgi and Manohar's naive dimensional analysis (NDA) |26U27| and 
naturalness, namely, that all appropriately defined, dimensionless couplings are of 
order unity. With this input, the nonlinear chiral Lagrangian can be organized in 
increasing powers of the fields and their derivatives. To each interaction term we 
assign an index 

n 

v = d+-+b, (2) 

where d is the number of derivatives, n is the number of nucleon fields, and b is the 
number of non-Goldstone boson fields in the interaction term. Derivatives on the 
nucleon fields are not counted in d because they will typically introduce powers of 
the nucleon mass M, which will not lead to small expansion parameters. 

It was shown in [19U281 I5] that for finite-density applications at and below nu- 
clear matter equilibrium density, one can truncate the effective Lagrangian to terms 
with v < 4. It was also argued that by making suitable definitions of the nucleon 
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and meson fields, it is possible to write the Lagrangian in a "canonical" form con- 
taining familiar noninteracting terms for all fields, Yukawa couplings between the 
nucleon and meson fields, and nonlinear meson interactions [2H]- See [51IU] for a more 
complete discussion. 



as 



If we keep terms with v < 4, the chirally invariant Lagrangian can be written 

2,3 



-Ceft = C-n + £4 + C-m 

= N{ij» [dp + iVf, + ig pPil + ig v V p ] + 34 7^75^ — M + g s 4>) N 



^LWv^N + ^-NN Tr (a^) + £ 4 



-a 2 J1tt(o v n + -(l + m^ + ^^)m 2 VV^ 
+ J { Co^Vf + (l + V„ f^) m 2 p Tr (p^) 

where the nucleon, pion, sigma, omega, and rho fields are denoted by N, 7T, 0, V^, 
and = I T-p p , respectively, V^j, = c^K, — c^V^, and a^ v = §[7^, 7"]. The trace 
"Tr" is in the 2x2 isospin space. The pion field enters through the combinations 

U = expiir-ir/U) , £ = exp(zr •ir/2f 7V ) , (4) 

^ = -^(^-^ f )=at, (5) 

v, = -~(^ + ^ + ) = ^' (6) 
*V = d p v u - d u v p + i[v p , v u ] = -i[a^ a v ] . (7) 



2 We use the conventions of [3]|SJ||5]- The pion-decay constant is fa 93MeV 
3 The two terms involving at coefiicier 
ically significant in See Fig.|SJ below 



3 The two terms involving at coefficients actually have v — 5, but they were found to be numer- 



Covariant Effective Field Theory 7 

The rho meson enters through the chirally covariant field tensor 

/V = D i*P" ~ D vP» + i 9p[Pn, Pv\ , (8) 
where the covariant derivative is defined by 

Dfj,p v = duPv + ifymPv] ■ (9) 



The antisymmetric combination of derivatives in p^ u implies that the timelike 
components Pq of the rho field have no conjugate momenta and are thus deter- 
mined by equations of constraint, as appropriate for a massive vector field with 
three dynamical degrees of freedom. The final term in (JSJ) has the usual form for a 
nonabelian vector field and enables the rho meson to couple to a conserved isovec- 
tor current [31120] • £4 contains tttt and 7rN interactions of order v = 4 that are 
not needed in this work. Numerically insignificant v = 4 terms proportional to 
4> 2 Tr (p M p M ) and V V V V Tr (p^p^) have been omitted. 

To exhibit the chiral invariance of £eft explicitly, we follow CCWZ |31 [l32[ l9j. 
A nonlinear realization of the chiral group S\J(2) L x SU(2)^ is defined such that for 
arbitrary global matrices L G S\J(2) L and R G SU(2)^, there is a mapping 

L®R: & P(1 ,N) -^(£',^0 . (10) 

Because of the parity operation V, which produces the transformation 

V , L^R, vr«(t,x) ^-vr a (t,-x) , £(t, x) — -x) , (11) 

the chiral mapping (|10j) can be written as jHJ 

£'(x) = L£(x)h\x) = h{x)i{x)R ] , (12) 
p'n(x) = h(x)p^{x)h ] {x) , (13) 
N'[x) = h(x)N(x) . (14) 

The second equality in (|12l) defines h(x) as a function of L, R, and the local pion 
fields: h(x) = h(L, R, tt(x)). It follows from l|12|) that h(x) is invariant under the 
parity operation pip, that is, 

h(x) G S\J(2) V , (15) 

where SU(2)^ is the unbroken vector subgroup of SU(2) L x SU(2) R . Note that the 
matrix h(x) becomes a constant only when L = R, so that h = L = R. Equa- 
tions pup and p4p then ensure that the rho and nucleon fields transform linearly 
under global SU(2)y, in accordance with their isospins. The isoscalar fields V^{x) 
and <f>(x) are chiral scalars and are unaffected by both chiral and isospin transfor- 
mations. 
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For discussing purely pionic interactions, it is convenient to use the matrix U(x) 
of @, since the transformation law ()12j) then implies 

U(x) — ► U'{x) = LU{x)R^ , (16) 

so that U(x) always transforms globally. Thus derivatives of U(x) transform the 
same way as U(x), and chirally invariant interactions involving pions alone can be 
constructed from products of U(x), W(x), and their derivatives. As is well known, 
these terms can be organized according to the number of derivatives, resulting in 
the Lagrangian of chiral perturbation theory [2*ll33|. We will return to this later 
when we discuss electroweak interactions with nuclei. 

For describing the interactions of pions with other particles, U (x) is not conve- 
nient, because other fields transform with the local functions h{x) of the unbroken 
isovector subgroup SU(2)y. It follows from the transformation laws given earlier 
that interaction terms that are invariant under local isospin rotations will be invari- 
ant under global transformations of the full group SU(2) L x SU(2)^. Thus, to form 
chirally invariant interactions involving pions and other fields, we need functions of 
the pion field that transform with h(x) only. 

The desired functions involving one derivative of the pion field are given in (|5|) 
and l[5|). The parity transformation implies that is an axial vector and i>„ is 
a polar vector. Moreover, under a chiral transformation, (|12j) implies 

a» — ► «V = ha^h) , (17) 
Vfj, — ► v'^ = hv^hJ — ihdfth* = hv^h) + i(d^h)h} . (18) 

Thus transforms homogeneously under the local SU(2)y group and can be in- 
terpreted as a covariant derivative of the pion-field matrix In contrast, the 
inhomogeneous transformation law for resembles that of a gauge field, so that 
allows us to construct chirally covariant derivatives of the other fields. For example, 
it is straightforward to verify that the covariant derivatives © and 

D„N = + iv^N (19) 

transform homogeneously with h(x) under the full group: 

(DpN)' = h(D^N) , (D^)' = h{D, lPv )h) . (20) 

The covariant tensor for the pion field is Vav [see (|Tjl]. which transforms homo- 
geneously with h, as does Plxv . This allows us to produce a chirally invariant pirir 
coupling through an interaction of the form Tr (p^ y v^ v ). 

Electromagnetic interactions can be included by adding a chirally noninvariant 
Lagrangian [we exhibit terms of 0(e) only] 
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£em = -\F^ -2eflA^(v^)-^-F^lTr{p^T Z ) + \V^ 



' d v F» v N 



2M 2 



(/3(°)+/3«T 3 ) 7/i ]iV , (21) 



where is the photon field and = e^-Aj/ — ^A^, is the usual Faraday tensor. 
The constants e, and (5^ in (|2~Tj) (where t = 0, 1 denotes the isospin), together 
with the tensor couplings f v and f p in (j3J , are sufficient to parametrize the empirical 
nucleon charge e, the anomalous moments A pn , and the charge and magnetic radii 

( r rms)i*2 a t l° w momentum transfer |Sj. The expansion can be extended to include 
higher derivatives of the photon field if greater accuracy is needed [21] • 

Similarly, the free parameter g& in the pion-nucleon interaction [see ©] allows 
us to normalize the one-body, axial-vector nuclear current so that the Goldberger- 
Treiman relation is satisfied at the tree level |H4] . 

To summarize the important points of the full Lagrangian £qhd = £eft + ^em 
[recall © and (HQ)]: 

• The noninteracting hadron terms take their standard canonical forms. 

• The generalized coordinates (fields) have been chosen so that the meson- 
nucleon couplings have a simple Yukawa form. 

• The pion-nucleon and pion-meson interactions enforce the nonlinear realiza- 
tion of chiral symmetry. 

• The nonlinearities involving chiral singlet fields are obviously invariant, and 
fitting their coefficients to data will implicitly include short-range dynamics 
from many-nucleon forces, fluctuations of the quantum vacuum, and hadron 
substructure. 



The nucleon electromagnetic (and weak) structure (gA, A, etc.) is included to 
the desired accuracy using a derivative expansion of the fields. 



3 Density Functional Theory 

The successes of QHD mean-field phenomenology are, at first, rather mysterious 
from the EFT perspective alone, since the Hartree approximation is just the finite- 
density counterpart of the Born approximation at zero density. The density func- 
tional theory (DFT) perspective explains the successes of mean-field approaches and 
provides a new context for EFT power counting. 
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We begin with a discussion of nonrelativistic DFT and generalize later to include 
relativity. The basic idea behind DFT is to compute the energy E of the many- 
fermion system (or, at finite temperature, the grand potential O) as a functional 
of the particle density [33] • DFT is therefore a successor to Thomas-Fermi theory 
[3011311 , which uses a crude energy functional, but eliminates the need to calculate 
the many-fermion wave function. 

The strategy behind DFT can be seen most easily by working in analogy to 
thermodynamics [38| . For a uniform system in a box of volume V at temperature T, 
one first computes the grand potential 0(/i, T, V), where is the chemical potential. 
It then follows that the number of particles N is determined by [Q 

"-<*>--(£!),,■ (22> 

According to Gibbs, thermodynamic equilibrium is defined by the condition 

(Sn)^ T ,v > ; (23) 

an assembly minimizes its grand potential at fixed fi, V, and T. Thus the convexity 
of f2 implies that N is a monotonically increasing function of fi, so the relation (j22j) 
can be inverted for n(N). Finally, one makes a Legendre transformation to the 
Helmholtz free energy 

F{N,T,V) = Sl(n(N),T,V)+n{N)N (24) 

to discuss systems with a fixed density n = N/V. 

For a self-bound, finite system, we replace the chemical potential with an ex- 
ternal, single-particle potential 4 Y^i v ( r i)- The grand potential is now a functional: 
fi[v(r);T), and a functional derivative with respect to v gives the particle density: 5 

n(r) = (n(r)) = ^ . (25) 

The convexity of O allows us to invert this relation (in principle) and to find v(r) 
as a (complicated) functional of n(r): 

v(r) = v[n{r)] . (26) 

(T is suppressed.) Thus there is a one-to-one relation between the external potential 
and the particle density. Needless to say, the possibility of this complicated inver- 
sion is a matter of great technical interest, and the reader is referred to [Six for a 
discussion. 



4 Here the chemical potential fi is absorbed into the definition of v, which defines the zero of 
energy. We suppress all spin and isospin dependence at this point. 
5 Higher variational derivatives yield various correlation functions. 
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Finally, we make a functional Legendre transformation to define the Hohenberg- 
Kohn free energy, which is a functional of n(r): 

F KK [n(r)] = Q [v[n(r)\ ] - J dr n{r)v{r) . (27) 

The variational derivative of this free-energy functional with respect to n now gives 

£-Fhk [n] 



8n(r) 



-v(r) , (28) 



where we have used (|25j). 

If we now restrict consideration to T = and v (r) = 0, the free energy be- 
comes simply the energy: -FhkM — ► ^hkM; an d the Hohenberg-Kohn theorem 
follows [9*ll39| : If the functional form of -EhkW 1 ")] is known exactly, the ground-state 
expectation value of any observable is a unique functional of the exact ground-state 
density. Moreover, it follows immediately from (|28|) that the exact ground-state 
density can be found by minimizing the energy functional. Although we have as- 
sumed here that the ground state is non-degenerate, this assumption can be easily 
relaxed |39| . 

Significant progress in solving these equations was made by Kohn and Sham 
|4Uj . who introduced a complete set of single-particle wave functions. The exact 
Hohenberg-Kohn free energy for an inhomogeneous (finite) many-body system in 
an external potential takes the form 

Fuk [n(r)\ = F ni [n(r)j + E int [n(r)] , (29) 

where the subscripts "ni" and "int" denote noninteracting and interacting, respec- 
tively. -F n i[n(r)] represents the kinetic energy contribution. The interaction energy 
£int[^(r)] is some functional of the density (and its derivatives); in the many-body 
problem, it contains a Hartree term, an exchange-correlation ("xc") contribution, 
etc. P: 

E int [n] = ^Hartree N + E xc [n] -\ . (30) 

Note that E xc is generally a nonlocal and nonanalytic functional of the density that 
contains both many-body and short-distance physics, including vacuum fluctuations 
and hadron substructure. 

Now consider the (nonrelativistic) Schrodinger equation in a potential v e s(r), 
which is designed to give the exact density n(r): 

(-^V 2 + Wcff (r))^(r) = e^(r), (31) 

N 

n(r) = ]ThMr)| 2 . (32) 
i = i 
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In this problem, v e g plays exactly the same role as the previous v, and thus the 
Hohenberg-Kohn equation (|28|) gives 

M^T = ~ Vcs{r) ■ (33) 

By taking the variational derivative of (|2H|) with respect to n(r), using (|2*%|) . and 

rearranging terms, we find 

v eS (r) = v(r) + . (34) 

Upon setting the external v(r) = 0, we obtain the effective potential to be used in 
the Kohn-Sham (KS) equations (|3*Tj) : 

5E int [n] 

«-M = -s^j- ■ (35) 

Thus, if we know the exact interacting energy as a functional of the density, 
we can reproduce the exact interacting density using a set of single-particle wave 
functions. Kohn calls these the "density-optimal" single-particle wave functions [22] 
(as opposed to Hartree-Fock wave functions, which are "total-energy optimal"). 

The generalization of DFT to relativistic systems is straightforward The 
energy -Ehk now becomes a functional of both the ground-state scalar density p s and 
the baryon four-current density B^. Extremization of the functional gives rise to 
variational equations that determine the ground-state densities p s and p B = Bq. 

These equations can again be simplified by following the Kohm-Sham approach. 
In the relativistic case, the complete set of single-particle wave functions allows us 
to recast the variational equations as Dirac equations for occupied orbitals. The 
single-particle Hamiltonian contains local, density-dependent, Lorentz scalar and 
vector potentials, even when the exact energy functional is used. Moreover, one can 
introduce auxiliary (scalar and vector) fields $(r) and W(r), which correspond to 
the local potentials and can therefore be identified as relativistic KS potentials. The 
auxiliary fields are determined by extremizing the energy functional, which gives rise 
to a Dirac single-particle Hamiltonian. The isoscalar part (for spherical nuclei) looks 
like 

ho = -ia ■ V + (3[M - $(r)] + W{r) , (36) 

where M is the nucleon mass and we define M* = M — <£. The resulting coupled 
differential equations resemble those in a relativistic MFT calculation jHHOj- Note 
that $ need not be simply proportional to the isoscalar, scalar field (p. In fact, 
$ could be proportional to (as in the Walecka model), or could be expressed as 
a sum of scalar and vector densities (as in relativistic point-coupling theories), or 
could be a nonlinear function of <f> (as in modern chiral EFT's). 
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The strength of the KS approach rests on the following theorem: 

The exact ground-state scalar and vector densities, energy, and chemi- 
cal potential for the fully interacting many-fermion system can be repro- 
duced by a collection of (quasi) fermions moving in appropriately defined, 
self-consistent, local, classical fields. 

The proof is again straightforward [SHI- Start with a collection of noninteracting 
fermions moving in an externally specified, local, one-body potential. The exact 
ground state for this system is known: just calculate the lowest-energy orbitals 
and fill them up. 6 Therefore, if one can find a suitable local, one-body potential 
based on an exact energy functional, the exact ground state of that system can 
be determined. But this potential is precisely the v e g(r) discussed above, which is 
obtained by differentiating the interaction parts of -Fhk with respect to the various 
densities. The resulting one-body potential will generally be density dependent and 
thus must be determined self-consistently. 

Several points are noteworthy. Since the single-particle basis constructed as 
described above is again "density optimal" , the exact scalar and vector densities are 
given by sums over the squares of the Dirac wave functions, with unit occupation 
probability. Moreover, since these densities are guaranteed to make the energy 
functional stationary [the external v(r) = 0], the exact ground-state energy is also 
obtained. The proof that the eigenvalue of the least-bound state is exactly the Fermi 
energy is given in [12]. Note, however, that aside from this association, the exact 
Kohn-Sham wave functions (and remaining eigenvalues) have no known, directly 
observable meaning. 

If one knows the exact functional form of the energy on the densities, one can 
describe the observables noted in the theorem exactly (and easily) in terms of the 
Kohn-Sham basis. Observables of this type are typically the ones calculated in rela- 
tivistic MFT's 43,44,8 . Moreover, it has been known for many years 3,45 that the 
mean-field contributions dominate the single-particle potentials at ordinary densi- 
ties. Thus, by parametrizing the energy functional in a mean-field (or "factorized" ) 
form, and by fitting the parameters to empirical bulk and single-particle nuclear 
data, one should obtain an excellent approximation to the exact energy functional 
in the relevant density regime. This is the key to the success of relativistic MFT cal- 
culations, as we will verify below, using the effective chiral Lagrangian constructed 
in Sect. El 

4 Naive Dimensional Analysis 

There is still an important point to be addressed: we must understand how to ex- 
tract the dimensional scales of each term in the Lagrangian, so that the remaining 

6 For simplicity, we assume that the least-bound orbital is completely filled, so the ground state 
is non-degenerate. 
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dimensionless constants can be checked for naturalness. A naive dimensional anal- 
ysis (NDA) for assigning a coefficient of the appropriate size to any term in the 
effective Lagrangian has been proposed by Georgi and Manohar |261I27| . This al- 
lows for a determination of both the dimensional scales associated with each term 
and for the inclusion of an overall dimensionless constant that can be used to adjust 
the strength. The basic idea of naturalness is that once the appropriate dimensional 
scales have been extracted, the overall dimensionless coefficients should all be of 
order unity. 

The NDA rules for a given term in the Lagrangian density are |27j : 

1. ) Include a factor of l/f n for each strongly interacting field. 

2. ) Assign an overall factor of / 2 A 2 . 

3. ) Multiply by factors of 1/A to achieve dimension (mass) 4 . 

4. ) Include appropriate counting factors (such as 1/n! for 4> n ). 

Here f n ~ 93 MeV is the pion-decay constant, and the breakdown scale A ~ 600 MeV 
is taken as the generic large-momentum-cutoff scale, which characterizes the mass 
scale of physics beyond the Goldstone bosons. 

As noted by Georgi |27], rule (1) simply assumes that the amplitude for produc- 
ing any strongly interacting particle is proportional to the amplitude f n for emitting 
a Goldstone boson. This is a reasonable assumption, since f n is the only natural 
scale. Thus, by dividing each field by f n , we should arrive at a factor of 0(1). 
Rule (2) can be understood as an overall normalization factor that arises from the 
standard way of writing the mass terms of non-Goldstone bosons. For example, one 
may write the mass term of a isoscalar, scalar field <j){x) as 

1 2/2 1 f2i2 m s ^ 

2™s</> = aT/2 > ( 37 ) 

where the scalar mass m s is treated as roughly the same size as A. By applying rule 
(1) and extracting the overall factor of / 2 A 2 , the remaining ratios are of 0(1). Since 
all terms will have the same overall scale factor / 2 A 2 , higher-order terms or terms 
with gradients of fields will be suppressed by powers of 1/A relative to the leading 
mass terms, as a result of "integrating out" physics above the scale A. (A simple 
example is the low-momentum expansion of a tree-level propagator for a heavy 
meson of mass mjj, which leads to terms with powers of d 2 /m 2 H .) It is precisely 
because of these 1/A suppression factors and dimensional analysis that one arrives 
at rule (3). The origin of the combinatorial factors in rule (4) is discussed in 

Applying these rules to a generic term in an effective Lagrangian involving the 
isoscalar fields and the nucleon field leads to (the generalization to include the pion, 
rho, and photon is straightforward) |46l l5] 
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where ip is a baryon field, V is any Dirac matrix, derivatives are denoted generically 
by d, and we have allowed for the possibility of chiral-symmetry-violating terms 
that contain the small parameter m n /A. The product of all the dimensional factors 
then sets the scale in terms of the pion-decay constant f n and the EFT breakdown 
scale A. The overall coupling constant C is dimensionless and of O(l) if naturalness 
holds. 

These scaling rules imply that a general potential for the scalar meson can be 
expanded as 

Vs = m ^ U + ¥M + 4!M^ + '-j ' (39) 

in agreement with the corresponding term in (J3J). Here we have included a factor of 
1/Jtt for each power of (f>; these factors are then eliminated in favor of g s ~ M/f n , 
which is basically the Goldberger-Treiman relation Factorial counting factors 
are also included, since the NDA rules are actually meant to apply to the tree-level 
scattering amplitude generated by the corresponding vertex |47l l5] . 

The naturalness hypothesis states that after the dimensional factors and appro- 
priate combinatorial factors are extracted, the overall dimensionless coefficients [C 
in (|38]> and K3, K4, ... in (|39]) ] should be of order unity. It should be clear, how- 
ever, that the preceding arguments are not a proof of naturalness, since unknown 
physical scales could generate unnaturally large coefficients. Moreover, some fitted 
constants may be unnaturally small, which often signals a symmetry of the theory 
that has not yet been identified. These caveats notwithstanding, without the nat- 
uralness hypothesis it is basically impossible to construct an effective Lagrangian 
with any predictive power. 7 Until one can derive the effective hadronic Lagrangian 
from QCD, the naturalness hypothesis must be checked by fitting to experimental 
data, as we will do in the following section. 

If truncations of the EFT Lagrangian determined by NDA and naturalness are 
valid, it should also be possible to determine the level of truncation that exhausts 
the information content of the input data. In other words, we should be able to 
verify that adding terms with higher values of v does not improve the fits to the 
empirical observables of interest. To this end, it is useful to consider the v = 5 
interactions 

1 Ss^ 3 2 ±2 1 1 5 r s < / )3 1 2t/ \ru 1 1 /- 9s4> 2fir iru\2 I *r\\ 

A = "ST 5 + 3! ??3 M^ , 2 mv ^ y + 4! Cl lF 5v0 ^ } ' (40) 

to check that their contributions are either negligible or can be absorbed into slight 
modifications of the parameters in the v < 4 Lagrangian £qhd- 

7 The assumption of renormalizability also leads to a finite number of parameters and well-defined 
predictions, but does so by imposing unnatural restrictions on the Lagrangian, namely, that many 
parameters are identically zero in the absence of relevant symmetry arguments. 
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5 Mean-Field Theory of Nuclear Structure 

The mean-field equations and energy resulting from £qhd in © and (|2*T|) can be 
derived straightforwardly. The interested reader is referred to [H], where the pro- 
cedures and results are discussed thoroughly. One important result is that due to 
the additional nonrenormalizable interactions in £em between the nucleon and the 
electromagnetic field, and also due to vector-meson dominance, the computed nu- 
clear charge density automatically contains the effects of nucleon structure, and it 
is unnecessary to introduce an ad hoc form factor. 

As we will show in Sect. the full MFT Lagrangian derived from £qhd has 
more than enough parameters to accurately describe the desired nuclear properties 
discussed in Sect.^ The more important question is whether the parameters fitted 
to nuclei are natural. In [SJ, the parameters were determined by calculating a set of 
observables {X^} for several spherical nuclei and by adjusting the parameters to 
minimize a generalized x 2 defined by \'2'A\ : 



where i runs over the set of nuclei, X runs over the set of observables, the subscript 

(i) 

"exp" indicates the experimental value of the observable, and W x are the relative 
weights. The weights were chosen to be the expected accuracy for the given ob- 
servable in a good fit. In practice, a reasonable range of weights was tested, and 
the qualitative conclusions discussed below were always reproduced. Some of the 
considerations relevant to the choice of weights are discussed in [HUH] . 

The relativistic mean-field equations are solved self-consistently for the closed- 
shell nuclei 16 0, 40 Ca, 48 Ca, 88 Sr, and 208 Pb, and also in the nuclear matter limit. 
The parameters are then fitted to empirical properties of the charge densities, the 
binding energies, and various splittings between energy levels near the Fermi surface 
using the figure of merit (x 2 ) in (|41|) . The full set of {X"} comprised 29 calculated 
and empirical values. When working at the highest order of truncation (essentially 
v = 4), the calculated results are very accurate, as we illustrate shortly, but they are 
too numerous to reproduce here [5ll24lllUj . Some of the fitted parameters at various 
levels of truncation are given in Table 




(41) 
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Table 1: Parameter sets from fits to finite nuclei, as described in the text. The 
parameters in the lower portion of the table are fitted to the (free) nucleon charge 
and magnetic form factors, and the proton charge e and the anomalous moments 
A p and A n are given their empirical values [Sj. 
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-0.38482 
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-0.01915 
-0.07120 


1.0332 
-0.10689 
-0.26545 


1.0660 
0.01181 
-0.18470 


1.0393 
0.02844 
-0.24992 


0.9619 
-0.09328 
-0.45964 
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To simplify the initial discussion, we restrict consideration to infinite nuclear 
matter. For symmetric matter (N = Z), the energy density through order v = 4 is 
given by [Hj 



£mft[$, W; Pb ] = Wp B + / F d 3 A: \/k 2 + M* 2 

( 27r r Jo 



gl \2 3! M 4! M 2 



1 



2 M 2 / 



miW z 



4<3 2 



Co^ 4 , (42) 



where fc F is the Fermi wavenumber, p B = 2kp/3ir 2 , and <1? = g s ^o = M — M* and 
W = g v Vo are the scaled fields defined earlier in terms of the scalar and vector mean 
fields (f)Q and Vq. For readers who are familiar with the corresponding result in the 
Walecka model (Eq. (3.53) in [H]), one sees that the MFT nuclear matter energy 
has been generalized to include additional nonlinearities that are not allowed in the 
renormalizable Walecka model. The fields <3? and W are determined by extremization 
of£. 

What causes the nuclear matter saturation and the relatively small binding en- 
ergy? Let's expand £mft/Pb from (|4*2*|) in powers of k F [S] and suppress the nonlinear 
meson terms for clarity: 



£mft/Pb 



M + 



3/c 2 



3/ctt 



+ 



10M 56M 3 48M 5 



gl 



m. 



_ PB 

2 M 



3A%j 



36^ 



10M 175M 3 



+ ... 



+ 



+ 



2m 2 



Pb 



g; 



9b Pb 



IM 



2m{ 



3/Ctt 



10M 



+ 



/ 5, 2 Pb V 


3k 2 


Vm 2 M / 


10M 



(43) 



The lowest-order Lorentz scalar and vector contributions (which are proportional 
to p B ) set the scale for the large mean fields $ and W. This scale is consistent 
with chiral QCD counting rules |8"lll0|. but these two terms cancel almost exactly in 
the binding energy, leading to an anomalously small remainder. However, they add 
constructively in the spin-orbit interaction, leading to appropriately large spin-orbit 
splittings in nuclei |4lfll4lfllTi] . 

It is important to notice the different behavior of the vector and scalar interaction 
terms in (|43j) . Whereas the vector interaction enters at only linear order in p B , 
the scalar interaction enters at all orders; moreover, the leading scalar term at 
every order in p B looks exactly the same, and they all add constructively. These 
terms are precisely what one gets by shifting the nucleon mass in the nonrelativistic 
kinetic energy term 3fc 2 /10M from M — » M* ~ M — g 2 p B /m 2 . These additional, 
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Figure 1: Nuclear matter energy /particle for two QHD parameter sets, one on the 
left (Gl) and one on the right (G2) of the error bars. The power of fields is b = j +£ 
for a term of the form Q^W^ (t is even). The boxes denote terms with j = 0, the 
circles denote terms with I = 0, and absolute values are shown. The crosses with 
error bars are estimates based on (|38[). with 1/2 < C < 2. The arrow indicates the 
total binding energy eo = 16.0 MeV. 



repulsive, velocity-dependent interactions reduce the strength of the lowest-order, 
attractive scalar contribution and are crucial for establishing the location of the 
equilibrium point of nuclear matter. Thus the different behavior of the vector and 
scalar interactions leads to large relativistic interaction effects in the nuclear matter 
energy density. In contrast, the relativistic corrections to the kinetic energy (the 
nonleading terms in the first pair of square brackets) are indeed small; this is not 
where the important "relativity" is. 

The critical question is whether the hierarchal organization of interaction terms 
is actually observed. This is illustrated in Fig. ^ where the nuclear matter en- 
ergy/particle is shown as a function of the power of the mean fields, which is called 
b in (J2J). (There are no gradient contributions in nuclear matter and (tt) = 0.) The 
crosses and error bars are estimates based on NDA and naturalness, that is, overall 
coefficients are of order unity. It is clear that each successive term in the hierarchy is 
reduced by roughly a factor of five, which implies a value of A ~ 600 MeV. Thus for 
any reasonable desired accuracy, the Lagrangian can be truncated at a low value of 
v. In fact, contributions to the energy /particle that are smaller than roughly 1 MeV 
are below the level of resolution and can be eliminated in favor of small adjustments 
of the remaining parameters. Derivative terms and other coupling terms will be 
discussed later. 

The quality of the fits to finite nuclei and the appropriate level of truncation 
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Figure 2: x 2 values for QHD parameter sets, as a function of the level of truncation. 
The power v = 2 corresponds to set Wl, v = 3 is for set CI, the u = 4 square 
is for set Ql, the v = 4 diamond is for set Q2, and the circle is for the full set of 
parameters Gl. The v = 5 results include the terms from £5 in (|40j) . 
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is illustrated in Fig. [2] 10j, where the figure of merit is plotted as a function of 
truncation order and of various combinations of terms retained in £qhd- The full 
calculations (•) retain all allowed terms at a given level of i/, while the other two 
choices keep only the indicated subset. There is clearly a great improvement in 
the fit (more than a factor of 35) in going from v = 2 to v = 4, but there is no 
further improvement in going to v = 5, using the extra interactions contained in 
(|40|) . Speaking chronologically, the v = 2 results show the level of accuracy obtained 
more than 20 years ago [13], while the v = 4 results were obtained seven years ago (Sj. 
Moreover, the li (p n only" results at v = 4 (■) show the state of the situation in the 
late 1980s, as discussed in |49l l7| . Recent work |10| shows that the full complement 
of parameters at order v = 4 is underdetermined, and that only six or seven are 
determined by this data set, which explains the success of the earlier MFT's with a 
restricted set of parameters {HI . We will review this analysis in the next section. 

As a further example of the quality of the fits, and for some additional historical 
perspective, Fig.|3]shows the percent deviation between the calculated and empirical 
binding energies for the five closed-shell nuclei listed earlier. The results labeled 
"Walecka" are from those labeled "QMC" are from recent quark-meson models 
of nuclear structure [230302]) and those labeled "QHD" follow from the present 
EFT for various parameter sets listed in Table ^ It is obvious that the modern 
EFT approach improves the quality of the fits by roughly two orders of magnitude 
and establishes a new standard of accuracy that must be attained for any modern 
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Figure 3: The deviations (in percent) between the calculated binding energy /nucleon 
and the empirical values for five doubly-magic nuclei. The different calculations and 
parameter sets are discussed in the text. For orientation, the Walecka- model results 
are underbound. 



approaches to nuclear structure to be considered viable. 

Finally, to illustrate the power of the EFT approach to nuclear many-body 
physics, we show two recent results of Huertas [S3]- The calculations use the MFT 
of £qhd with parameters fitted to closed-shell nuclei along the "valley of stability" , 
namely, set G2 in Tabled The resulting self-consistent, relativistic Kohn-Sham and 
meson field equations are solved to calculate the properties of Sn isotopes out to 
doubly-magic values of N and Z far from stability. These results are therefore true 
predictions of the EFT, since no adjustments were made to previously determined 
parameters. Figure^shows the predicted binding energies of the even-even isotopes 
from 1 5(]Sn to 1 |oSn compared with measured experimental results. The theoretical 
values are accurate to better than 1% throughout. (Similar accuracy is obtained for 
parameter set Gl; see Fig. 4 in [53*].) 

Figure shows the predicted ground-state energies, spins, and parities of the 
neighboring single-particle and single-hole nuclei 1 |iSb and 1 |gln relative to foSn. 
The energy differences are just the chemical potentials, which should be accurately 
reproduced, according to the discussion of Kohn-Sham theory in Sect. 01 The 
excellent agreement provides compelling evidence that QHD is indeed an EFT for 
low-energy QCD that can be used to describe nuclear many-body physics. This 
recent work has been extended to semi-magic nuclei with N = 28,50,82,126 and 
Z = 28,50,82 in 

How can we understand the excellent accuracy of the preceding MFT results? 
As discussed in Sect. [3 the exact energy functional has kinetic-energy and Hartree 
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Figure 4: Binding energy of even-even 5 ^Sn isotopes calculated using the MFT of 
£qhd with parameter set G2 from Table ^ „53 j . 
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parts (which are combined in a relativistic formulation) plus an exchange-correlation 
functional, which is generally a nonlocal, nonanalytic functional of the densities that 
contains all the other many-body, relativistic, and short-range effects |49[l55j . The 
basic idea behind the relativistic MFT (RMFT) is to approximate the functional us- 
ing an expansion in classical meson fields (or nucleon densities) and their derivatives, 
based on the observation that the ratios of these quantities to the nucleon mass are 
small, at least up to moderate density. 8 The parameters introduced in the expansion 
are fitted to experiment, and if we have a systematic way to truncate the expansion, 
the framework is predictive. Moreover, if the RMFT energy functional is sufficiently 
general, it will automatically incorporate effects beyond the Hartree approximation, 
such as those due to short-range physics and many-body correlations. 

But why should we expect an approximate, mean-field energy functional to work 
so well? We observe that while the mean scalar and vector potentials $ and W are 
small compared to the nucleon mass, they are large on nuclear energy scales |56II19| . 
Moreover, as illustrated in Dirac-Brueckner-Hartree-Fock (DBHF) calculations |45[ 
157118] . the scalar and vector potentials (or self-energies) are nearly state independent 
and are almost equal to those obtained in the Hartree approximation. Thus the 
Hartree contributions to the energy functional should dominate, and an expansion 
of the exchange-correlation functional in terms of mean fields should be reasonable. 
This "Hartree dominance" also implies that it should be a good approximation 
to associate the single-particle Dirac eigenvalues with the empirical nuclear energy 
levels, at least for states near the Fermi surface |35|. 9 

We also observe that the nuclear properties of interest discussed in Sect. E 
include: 1) nuclear shape properties, such as charge radii and charge densities, 
2) nuclear binding-energy systematics, and 3) single-particle properties such as level 
spacings and orderings, which reflect spin-orbit splittings and shell structure. Since 
the Kohn-Sham approach is formulated to reproduce exactly the ground-state en- 
ergy and density, and the Hartree contributions are expected to dominate the Dirac 
single-particle potentials, these observables are indeed the ones for which meaningful 
comparisons with experiment should be possible. 

As discussed above, an RMFT energy functional of the form in (|42|). extended 
to include low-order derivatives of the meson fields, successfully reproduces these 
nuclear observables with parameters of natural size (see Table This justifies a 
truncation of the energy functional at the first few powers of the fields and their 
derivatives, as is evident from Fig.0 Moreover, the full complement of parameters is 
underdetermined, so keeping only a subset does not preclude a realistic fit to nuclei. 

8 Since the meson fields are roughly proportional to the nuclear density, and since the spatial 
variations in nuclei are determined by the momentum distributions of the valence-nucleon wave 
functions, this organizational scheme is essentially an expansion in k F /M, for fc F corresponding to 
ordinary nuclear densities. Here the nucleon mass M ~ A is a generic large mass scale characterizing 
physics beyond the Goldstone bosons. 

9 One expects the KS spin-orbit splittings to be more accurate than the absolute energy eigen- 
values. 
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Both the early RMFT calculations mentioned above and the newer calculations 
based on chiral EFT should be interpreted within the context of this Kohn-Sham 
approach to DFT. 

6 Analysis of Mean-Field-Theory Parameters 

Although we could use meson-nucleon EFT's as done historically, the analysis is 
more transparent with point-coupling theories, which contain only nucleon fields in 
a local Lagrangian. Because of the freedom to perform field redefinitions, a general 
point-coupling Lagrangian is equivalent to a general meson-nucleon Lagrangian [SI 

C2H23- 

An energy functional of nucleon densities can be constructed by starting with a 
general point-coupling effective Lagrangian, consistent with the symmetries of QCD, 
and by evaluating the corresponding one-loop energy functional. As discussed above, 
this approach approximates a general DFT functional that incorporates many-body 
effects beyond the Hartree level when the parameters are determined from finite- 
density data. 

To arrive at a suitable truncation scheme, we again rely on NDA and naturalness. 
The two relevant mass scales are f n and A, and for closed-shell nuclei, the energy 
functional is an expansion in powers of the nucleon scalar, vector, isovector- vector, 
tensor, and isovector-tensor densities, which are defined as p s = (ipip), Pb = (ip^ip), 
P3 = ^{^tsi/j), s i = (tpo' 0t ip), and = ^(ip<J 0i Tsil'), respectively, where ip is the 
nucleon field. 

We can then define scaled densities and their derivatives as [see (J5S)) ] 

^ s " if/V ' = ' etC ' (44) 

NDA also provides numerical estimates for the scaled densities that will allow us to 
estimate terms in the energy functional. For example, each additional power of p s 
is accompanied by a factor of /^A. The ratios of scalar and vector densities to this 
factor at nuclear matter equilibrium density are between 1/4 and 1/7 j^Hj, which 
serves as an expansion parameter. Similarly, one can anticipate good convergence 
for gradients of the densities, since the relevant scale for derivatives in finite nuclei 
is the nuclear surface thickness <t, and so the dimensionless expansion parameter 
is 1/Acr < 1/5. The expansion is useful because the coefficients have been shown 
empirically to be natural, that is, of order unity [§ll24| . 

The energy functional is dominated by the isoscalar terms, but we also include 
the isovector and tensor terms for completeness. To nominal order v = 4 (all 
densities and gradients count as v = 1, except the three- vector tensor densities s, 
which have v k, 2; see also footnote |3] on p. El), we obtain 
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where the notation of [21] is used. The sum runs over occupied nucleon states. 

In Fig. the small symbols show the NDA estimates (with associated error 
bars) for the various energy contributions in 16 O and 208 Pb. The magnitudes of 
energy contributions in 1)451) from two representative RMF point-coupling models 
(i.e., two different parameter sets) are shown as larger unfilled symbols (one model 
on each side of the error bars). These models provide very accurate predictions 
of bulk nuclear properties [23]. The energy contributions are determined for each 
nucleus by making multiple runs while varying each parameter slightly around its 
optimized value, which enables us to deduce the logarithmic derivative with respect 
to each parameter. The filled symbols denote the sum of the values for each power 
of the density. The binding energy/nucleon in equilibrium nuclear matter is denoted 
by e . 

The two representative point-coupling models validate the isoscalar estimates 
(small open squares), and the resulting hierarchy of isoscalar contributions is quite 
clear. How far down in the hierarchy can we reliably determine contributions and 
their associated parameters? In Fig. [2] the impact of different truncations of RMF 
meson-nucleon models is shown by plotting the figure of merit against the maximum 
power of fields. We have also performed this test with point-coupling models, with 
similar conclusions. The "full" models (which include all nonredundant terms at a 
given order) show that one needs to go to the fourth power of the fields to get the 
best fits, but going further yields no improvement. Analogous behavior is found for 
RMF point-coupling models with powers of densities replacing powers of fields [51] . 
Thus contributions to the energy /particle at the level of roughly lMeV are at the 
limit of resolution. Fifth-order isoscalar contributions to the energy/particle, which 
are predicted to be less than 1 MeV, are simply not determined by the optimization. 

The variation in coefficient values in (|45|) provides a measure of how well the 
parameters are actually determined by the data. Figure shows the seven coeffi- 
cients of isoscalar non-gradient terms from four point-coupling models. Note that 
all coefficients are natural, i.e., of order unity. However, the spread in coefficient 
values is significant and does not correspond well to the power-counting order. We 
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Figure 6: Contributions to the energy /particle in 16 and 208 Pb determined by 
logarithmic derivatives with respect to the model parameters (see text) for two 
RMF point-coupling models 24 . Absolute values are shown. The filled symbols 
are net values. The small symbols indicate estimates based on NDA ^U^, with the 
error bars corresponding to natural coefficients from 1/2 to 2. 
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Figure 7: Coefficients for four accurately fit RMF point-coupling models from |24j . 
Each model is represented by a different shape, and the shading shows the type of 
term (scalar, vector, or mixed). 



Table 2: Improved coefficients for point-coupling RMF models (Table VI of [24J. 





linear 


density 


deduced 


coefficient 


combination 


scaling 


value 


n L 


«2 + C2 


P+ 


-0.51 ±0.01 


n 3 


K 3 + 7/1 


Pi 


+1.3 ±0.1 


n 2 


K 2 ~ C2 


P- 


-2.0 ±0.4 


n 5 


k 4 + C4 + m 


P\ 


-2.4 ±0.7 


n 4 


^3 - 7/1/3 


P+P- 


+0.2 ±1.0 


n 6 


K4 - C4 


p\p- 


-2.6 ±0.8 


n 7 


k 4 + C4 - 2r/ 2 /3 


P+P- 





conclude that different linear combinations of the coefficients must be considered to 
draw reliable conclusions about how many are determined by the data. 

Can we find a more systematic power counting scheme? The similar size of the 
scalar density p s and the vector density p B suggests that we count instead powers of 
p + = (p s + /9 B )/2 and p^ = (p s — /9 B )/2. The corresponding "improved" coefficients 
are listed in Tabled [see (|45|)]. The spread in these coefficients for four RMF point- 
coupling models from [21] is shown in Fig. |H1 The terms are organized according to 
the powers of p + and p_, with /?_ scaling as /?+ 3 . 

The leading orders are very well determined, with a systematic increase in un- 
certainty. Even the sign is undetermined for the parameter Q.^, which is shown 
with unfilled symbols, but the next parameter (f^) appears to be reasonably well 
determined. Higher-order terms are not determined by the optimizations. Deduced 
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Figure 8: Improved coefficients for the same four models as in Fig. [7J The "order" 
is determined by counting powers of /?+ and p- . 

values and uncertainties based on this sample of models are given in Table We 
see that of the seven isoscalar non-gradient parameters in ([45 [I. four linear combi- 
nations are clearly determined by bulk nuclear observables, with probably a fifth 
combination as well. 

The isovector terms appear only on the graphs for 208 Pb. The factor (N—Z)/2B, 
which is only 10% even for Pb, severely limits the sensitivity to isovector terms 
(especially since the factor must appear with even powers). The magnitude of the 
leading isovector term (oc £,2f>i ) is comparable to the fourth-order isoscalar term, 
which is at the limit of what can be determined reliably from fitting the binding 
energy. We conclude that only one isovector parameter is determined by the bulk 
observables. For a further discussion of the isovector terms, see |59) . 

The energy estimates for isoscalar, tensor terms imply that only one parameter, 
at best, can be determined. Higher-spin terms, which will require more gradients 
and have smaller average densities, are not at all constrained. The tensor terms are 
interesting because a fraction of the spin-orbit force can be generated by including 
an isoscalar, tensor coupling of the vector field to the nucleon Jl]. Nevertheless, 
the spin-orbit potential arises predominantly from the large scalar and vector fields; 
attributing more than one-third of the potential to the tensor coupling produces 
unrealistic surface systematics 60 . 

Finally, gradient terms follow the same pattern in the energy: the leading term 
is barely above the limit of resolution. In fact, there are two isoscalar gradient terms 
at leading order (scalar and vector), but only their sum is well determined. The 
sum of the subleading-order contributions almost vanishes. We conclude that only 
one gradient parameter can be determined 59 . 

The handful of parameters that are well determined by the usual bulk nuclear 
observables (binding energies, charge density distributions, and spin-orbit splittings 
in doubly magic nuclei) can be associated with an equal number of nuclear properties 
and general features of RMFT's. In particular, 
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1. ) Two isoscalar non-gradient parameters are very well determined. These cor- 

respond to the highly constrained values for the equilibrium density (k F ~ 
1.30 ± 0.01 fm" 1 ) and binding energy (16.0 ±0.1 MeV) of nuclear matter. 

2. ) An additional isoscalar constraint is that M* ~ 0.61 ± 0.03, if the isoscalar, 

tensor term is set to zero. This range ensures an accurate reproduction of 
spin-orbit splittings in finite nuclei. Small increases in M* without changing 
the splittings can be accomplished by including an isoscalar, tensor term; an 
analysis using a simple local-density approximation is discussed in |14j . 

3. ) A fourth isoscalar constraint comes from the nuclear matter compressibility. 

The constraint is much weaker, in the range of K ~ 250 ± 50 MeV. 

4. ) The possibility of a fifth isoscalar constraint has been considered by Gmuca 

|61j , who argued that separate scalar and vector fourth-order terms were nec- 
essary to tune the density dependence of the scalar and vector parts of the 
baryon self-energy. This would correspond to constraining S7g from Table El 
Moreover, some form of isoscalar nonlinear vector interaction is needed to 
soften the high-density equation of state to be consistent with observed neu- 
tron star masses |55j . 

5. ) Since only one isoscalar gradient parameter is determined, it is not useful to 

allow the scalar and vector masses (or their equivalents in a point-coupling 
theory) to vary independently. Thus it is convenient to fix the vector mass 
at a natural size, such as the experimental mass for the u. A scalar mass of 
500 ± 20 MeV is then required. 

6. ) The one isovector parameter can be fixed by the surface-corrected volume 

symmetry energy |62) . which falls in the range 34 ± 4 MeV |43| . Since no 
isovector gradient is determined, setting the isovector, vector meson mass to 
the experimental p meson mass is adequate. 

To ensure a reasonable (if not optimal) description of finite nuclei, it is sufficient 
to reproduce the nuclear matter properties given above, but note that all calculated 
properties must be accurate, and the resulting parameters must be natural. 10 One 
cannot justify the underlying physics of a model if it reproduces only a subset of 
the nuclear calibration data. 



10 A further caution is that there are many correlations among these properties, so that the 
allowed ranges should not be considered to be independent. 
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7 Weak Nuclear Currents 

A desirable theory of nuclear currents should satisfy the following three criteria: 

• It should use the same degrees of freedom to describe the currents and the 
strong-interaction dynamics. 

• It should respect the same internal symmetries, both discrete and continuous, 
as the underlying theory of QCD. 

• Its parameters can be calibrated using strong-interaction phenomena, like 7rN 
scattering and the properties of finite nuclei. This is especially important in 
EFT's, as they contain all (non-redundant) interaction terms that are consis- 
tent with the underlying symmetries |63l l§] . 

The QHD framework described so far embodies these three desirable features. 

The weak currents arise from Noether's theorem applied directly to £qhd and 
contain the pion field to all orders. The leading-order (in v) vector and axial-vector 
currents are given by (a is the isospin index) 



& Tr {r a (udPtf + £/ f } + ~ 



1 ^ 



1 



N 



N 



N 



1 



N . 



(46) 



(47) 



As discussed in [Mj, the Feynman rules for the lowest-order 7rN vertices can be 
determined from (J3J), and in the presence of an external axial- vector source, one 
can compute the one- and two-nucleon amplitudes for the axial- vector current. The 
amplitude for axial-current pion production on a single nucleon is represented by 
the diagrams of Fig. El and can be written as 
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(48) 



where is the outgoing four-momentum of the emitted pion, and p^ and p ,fl are 
the initial and final nucleon four-momenta, respectively. 



Covariant Effective Field Theory 



31 



ua ® 



|j,a 



b ® 



/S + 

k 

n 



Figure 9: The amplitude for axial-current pion production on a single nucleon 
[Eq. (|48|)]. Note that momentum k^ is extracted by the external source (which 
is denoted by <g>), so that p 1 ' = p' M + + q^. 
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Figure 10: Two-nucleon contributions to the axial current. Here the vertex M(tt) 
[Eq. (|48() ] represents all possible ways that the external source can couple to the 
left-hand nucleon line (to the order we are working), as illustrated in Fig. |H1 There 
are also contributions from the right-hand nucleon line, as well as exchange terms. 



The two-nucleon amplitude for the axial current, to this order in v, is given 
in Fig. I1UI It is easy to verify that this amplitude satisfies PCAC (for on-shell 
nucleons). 

One can also determine the vector and axial-vector currents arising from the 
v = 3 and v = 4 terms in the QHD Lagrangian These additional contributions 
take the form 
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Figure 11: Two- nucleoli, axial-current amplitudes originating from the v = 3 pion- 
nucleon terms in the effective Lagrangian. There is one pair of diagrams proportional 
to K n and another pair proportional to (5^. 
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The terms proportional to n n and (3^ produce the two-nucleon amplitudes shown in 
Fig. 1111 The corresponding amplitudes involving rho meson exchange are illustrated 
in Fig. 1121 and there are also amplitudes with isoscalar scalar and vector meson 
exchange that resemble the first two diagrams in this figure. Analytical expressions 
for all of these two-nucleon amplitudes are given in [HI] . 

The scattering amplitudes constructed from the currents listed above satisfy 
CVC, PCAC (when m n ^ 0), and the Goldberger-Treiman relation (with ^ 1) 
automatically. Moreover, the chiral charges Q a and Q§ derived from these currents 
satisfy the familiar chiral charge algebra to all orders in the pion field |34j . In |64| , 
these currents are used to study beta decay in 131,133 Sn. 



8 Summary 

In this work, I discussed recent progress in Lorentz-covariant quantum field theories 
of the nuclear many-body problem, often called quantum hadrodynamics (QHD). 
QHD is a local, nonrenormalizable, effective Lagrangian field theory with baryons 
and mesons as the generalized coordinates (fields). An effective Lagrangian con- 
sists of known long-range interactions constrained by symmetries and a complete, 
non-redundant set of short-range interactions. By simply looking at the spectra 
of massive nuclei, it is obvious that some relativistic effects must be important in 
nuclei; thus, it is most convenient to use a Lorentz-covariant theory. 

The effective field theory studied here contains nucleons, pions, isoscalar scalar 
(a) and vector (u>) fields, and isovector vector (p) fields. The heavy mesons are 
introduced as collective, effective degrees of freedom to simplify the description 
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Figure 12: Two-nucleon, axial-current amplitudes containing rho meson exchange. 
The notation in the first two diagrams implies that the external source can couple 
to the nucleon line in two ways, as in Fig. 1111 



of the medium- and short-range nucleon-nucleon interaction and to conveniently 
parametrize ground-state expectation values of nucleon bilinears, which are impor- 
tant for the description of bulk nuclear properties. The QHD theory exhibits a 
nonlinear realization of spontaneously broken S\J(2) L x SU(2)^ chiral symmetry 
and has three desirable features: it uses the same degrees of freedom to describe the 
nuclear currents and the strong-interaction dynamics, it respects the symmetries of 
the underlying theory of QCD, and its parameters can be calibrated using strong- 
interaction phenomena. Moreover, the electromagnetic structure of the nucleon can 
be included straightforwardly in a derivative expansion of the fields. 

Although the QHD Lagrangian in principle contains an infinite number of terms, 
naive dimensional analysis and naturalness allow one to identify suitable expansion 
parameters and to estimate the sizes of various terms in the Lagrangian. Thus, for 
any desired accuracy, the Lagrangian can be truncated to a finite number of terms. 
In particular, for normal nuclear systems, it is possible to expand the QHD effective 
Lagrangian systematically in powers of the meson fields (and their derivatives) and 
to truncate the expansion reliably after the first few orders. 

Using density functional theory, I showed that the mean-field approximation 
produces an energy functional whose parameters can be determined by fitting bulk 
and single-particle properties of nuclei. The framework of Kohn-Sham theory allows 
the ground state to be constructed from (quasi)particle orbitals with unit occupation 
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number. Since the mean-field energy functional is a good approximation to the 
exact energy functional over the relevant range of density, it is possible to reproduce 
nuclear densities, binding energies, and single-particle spectra near the Fermi surface 
very accurately. Because the parameters are fitted to nuclear properties, the energy 
functional implicitly contains effects that go beyond a simple Hartree approximation, 
such as short-range physics, hadron substructure, and many-body correlations. 

The numerical parameters of QHD were studied using an effective, point-coupling 
Lagrangian that contains nucleon fields only. Because of the freedom to redefine the 
fields (or coordinates), a general point-coupling theory is equivalent to a general 
baryon-meson theory. By examining the contributions to the energy /nucleon in 
doubly magic nuclei, it was found that only a small number of parameters (roughly 
seven) can be calibrated by the nuclear data input. New ways to calibrate additional 
parameters will play an important role in the construction of the next generation 
of QHD Lagrangians. Finally, the weak vector and axial-vector currents, and the 
corresponding two-nucleon, axial-current amplitudes, were discussed in the QHD 
framework. 

Nuclear physics is the study of strongly interacting hadronic matter, and the only 
consistent theoretical framework we have for describing such a relativistic, interact- 
ing, quantum-mechanical, many-body system is relativistic quantum field theory 
based on a local Lagrangian density. Although QCD of quarks and gluons provides 
the basic underlying theory, Lagrangians comprised of hadronic degrees of freedom 
(QHD) provide the most efficient description of the physics in the strong-coupling 
domain. In the modern effective field theory perspective of QCD, one incorporates 
only the underlying symmetries of QCD into the QHD Lagrangian. By interpreting 
the mean-field approximation in the context of density functional theory, one can 
understand the numerous successes of the QHD description of nuclear properties. 
Nevertheless, finding an efficient, tractable, nonperturbative way to match the QCD 
Lagrangian to the long-range, strong-coupling, effective field theory of QHD is still 
a major goal for the future. 
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